Flow curves of dense colloidal dispersions: schematic model analysis of the 
shear-dependent viscosity near the colloidal glass transition 
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A recently proposed schematic model for the non-linear rheology of dense colloidal dispersions 
is compared to flow curves measured in suspensions that consist of thermosensitive particles. The 
volume fraction of this purely repulsive model system can be adjusted by changing temperature. 
Hence, high volume fractions (<j> < 0.63) can be achieved in a reproducible manner. The quantitative 
analysis of the flow curves suggests that the theoretical approach captures the increase of the low 
shear viscosity with increasing density, the shear thinning for increasing shear rate, and the yielding 
of a soft glassy solid. Variations of the high shear viscosity can be traced back to hydrodynamic 
interactions which are not contained in the present approach but can be incorporated into the data 
analysis by an appropriate rescaling. 
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I. INTRODUCTION 



The flow behavior of concentrated suspensions under 
steady shear is a classical subject of colloid physics^. A 
large number of experimental studies conducted mostly 
on hard spheres has established the basic facts: If the 
concentration of particles is not too high, a first Newto- 
nian region is observed if the shear rate 7 is small. Here, 
the solution viscosity 770 measured in this first Newtonian 
regime can be significantly larger than rj s the one of the 
pure solvent. At higher shear rates, the perturbation of 
the microstructure of the suspension by the convective 
forces can no longer be restored by the Brownian motion 
of the particles. Hence, significant shear thinning will 
result in which the reduced viscosity r]jr\ s is more and 
more lowered until (often) a second Newtonian region is 
reached. In this region, r]jr\ s is generally considered to 
be dominated by the hydrodynamic interactions between 
the particles. Highly concentrated suspensions behave as 
weak amorphous solids and elastically withstand small 
but finite stresses^. 

The experimental results obtained so far have demon- 
strated that the deviation from the equilibrium struc- 
ture can be gauged in terms of the bare Peclet number 
Pco = a 2 7/Z?o, where a denotes the particle radius and 
Do the diffusion coefficient at infinite dilution 1 . How- 
ever, non-Newtonian flow behavior is observed already at 
rather small Peo and the disturbance of the microstruc- 
ture sets in at a shear rate defining a second characteristic 
number, the Weissenberg number or dressed Peclet num- 
ber Pe=7T, which is connected to the structural relax- 
ation time t. Shear thinning may be considered to arise 
for Pe > 1 . In dispersions able to order under shear, the 
viscosity is observed to decrease. Yet, shear-thinning in 
concentrated suspensions of colloidal particles at low 7 is 
not necessarily related to the onset of crystallization or 
other effects solely occuring at very high shearing fields. 
The structure remains amorphous during the application 



of shear rates that already lead to a marked decrease of 
the shear viscosity^^. 

The enormous raise of the zero-shear viscosity 770 with 
increasing volume fraction has been a long-standing prob- 
lem in the fieldi. Earlier theoretical approaches^ have 
assigned this increase of 770 /n s to the onset of the struc- 
tural arrest if the system is approaching the volume frac- 
tion of random close packing located at ca. cf> = 0.63. 
Hence, rjQ/rj s is predicted to diverge at this limit. How- 
ever, Meeker et al. in 1997— carefully re-analyzed all ex- 
perimental data of 770 / rj s available at that time and con- 
cluded that the strong raise of the zero-shear viscosity 
is related to the glass transition in suspension occurring 
at the volume fraction <j) g « 0.58. Other comparisons 
strongly supported this view^, while recent viscosity mea- 
surements remained inconclusive 10 ' 11 . 

Nagele and coworkers worked out a theoretical 
approach^^iii that could explain the increase of 770 /rj s on 
the base of the mode coupling theory (MCT) of Gotze 
and coworkersi^i^. Hence, the marked slowing down 
of the mobility of concentrated suspensions can directly 
be traced back to the caging of a given sphere by its 
surrounding neighbors. The quantitative description of 
the dynamics of quiescent suspensions in terms of the 
MCT has met with gratifying success when confronted 
with experimental data obtained through dynamic light 
scattering 1 ?^ 1 ?! 1 ?! 20 ! 21 . Indeed, MCT was shown to ex- 
plain the structural arrest of concentrated suspensions 
and describe quantitatively the dynamics that stretch out 
over many orders of magnitude. Moreover, in an impor- 
tant paper Mason and Weita 2 ^ could demonstrate that 
MCT leads to a full explanation of the linear viscoelastic 
behavior of hard sphere suspensions near the glass tran- 
sition. 

Recently, a theoretical model for the shear-thinning of 
concentrated suspensions was presented 2 ^!. It is based 
on MCT and gives a full description of the reduced viscos- 
ity 77/775 as the function of the shear rate 7. A comparison 
with recent simulations has demonstrated that this the- 
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ory captures all the salient points of the flow behavior of 
glassy systems^ 5 -. In particular, theory predicts a finite 
yield stress beyond the glass point. It vanishes discon- 
tinuously when going below the glass transition, where 
a first Newtonian plateau appears, which is followed by 
strong shear thinning. 

In this paper we present the first comprehensive com- 
parison with experimental data obtained from a model 
system. The paper is organized as follows: In section II 
we review briefly the central concepts of our approach. 
Section III then gives the quantitative comparison with 
recent experimental data obtained on a model system2&. 
A final section will conclude this paper. 



II. FLOW CURVES AND THE COLLOIDAL 
GLASS TRANSITION 

A. Loss of structural memory caused by shear 
advection 

As mentioned above, the marked shear thinning in 
dense dispersions, that is, the speeding up of structural 
relaxation through shear is not necessarily related to 
shear ordering. This was shown by experiments^ as well 
as by Brownian dynamics simulations 27 . In Refsi2i24 we 
argued that the speed up of decorrelation brought about 
by shear advection combined with local Brownian motion 
lies at the origin of shear thinning in dense dispersions. In 
this contribution we work out the involved loss of struc- 
tural memory caused by shearing in a schematic model 
that captures the universal aspects of the full microscopic 
approach of Ref^. 

The approach of Refsi2i24, connecting the nonlinear rhe- 
ology of dense dispersions to the glass transition, predicts 
a transition from a shear-thinning fluid to a yielding 
solid. Even though small shear rates are considered and 
the (bare) Peclet number Peo = 7a 2 /D is negligible, the 
final relaxation of transient density fluctuations or of the 
transient stress moduli is strongly accelerated by shear 
whenever Pe= jr is not negligible. In fluid states, where 
t is large at 7 = 0, shear advection speeds up the decay of 
structural correlations. For states which would be solid 
without shear, and where the shear modulus G(t) would 
arrest at a (finite) elastic constant at long times, enforc- 
ing stationary shear leads to a finite relaxation tim o 23 i 24 
which is of the order of l^l" 1 . Hence, the glassy state of 
the suspension is shear-melted. The suspensions yields 
as stress fluctuations decay to zero with rate set by the 
external drive. 



B. Universal aspects 

The flow curves a versus 7 exhibit qualitative aspects 
that are solely determined by the nature of the transition. 
With the separation parameter e denoting the (relative) 
distance from the transition, and to the time scale ob- 



tained by matching onto microscopic short-time motion, 
the following behaviors of the steady state shear stress a 
in the 'structural window' have been established 25 

r 7*0 (-e)^G C oo e«-i7; | 1/7 

a = a(jt ,e) -» I a+ (1 + c 3 \^t \ m ) \e\ « \jt \& , 
{ <t+ (1 + c 4% /e) e>|7io|^ 

(1) 

where the appearing constants are positive material- 
dependent parameters and the exponents 7, a, and m 
are non-universal numbers that are uniquely determined 
by the quiescent static structure factor^ 2 *!. The first 
line of Eq. , which describes the divergence of the vis- 
cosity, is familiar from classical MCT, and is discussed 
in e.g. Refi 2 ^". The 'structural window', here, is de- 
fined as the double regime |e| <C 1 and |7i | -C 1, 
where the slowing-down of the structural dynamics dom- 
inates the steady state stress. A 'dynamic yield stress' 
<7 + (e) = (7(7 — > 0, e > 0) is obtained in the glass be- 
cause a finite stress has to be overcome in order to force 
the glass to yield even for vanishingly small shear rate. 
The given asymptotes are only the leading orders for 
e — ► and 7*0 — > 0, while corrections can be obtained 
systematically 2 ^^. A model calculation will be shown 
further below. 



C. Schematic models 

The universal phenomena summarized in Eq. l[TJl 
exist in any model that exhibits the bifurcation sce- 
nario from yielding solid to shear thinning fluid. The 
central feature of the equations of motion is that they 
contain the competition of two effects: i) a non-linear 
memory effect increases with increasing particle inter- 
actions ('collisions' or 'cage effect') which leads to a 
non-ergodicity transition in the absence of shear, and 
ii), memory effects vanish with time because of shear- 
induced decorrelation. Both effects can be captured 
in the simpler 'schematic' models also. Note that the 
models can be set-up so that they obey similar stability 
equations as the microscopic approach. Thus, the corre- 
sponding asymptotic results summarized in Eq. Q hold. 

The well studied and comparatively simple schematic 
F^-model considers one normalized correlator $(t), 
which obeys a generalized relaxation equation 2 *: 

$(*) + T j$(t) + Jdt' m(t - t') $(i') J = . (2) 

Without memory effects, m = 0, the correlator relaxes 
exponentially, <fr(t) — exp— Tt, but with m =^ 0, retarda- 
tion effects set in after a short-time variation (still given 
by the initial decay rate T, viz. $(t — ► 0) = 1 — Tt + . . .). 
The correlator <I>(t) is taken to model the normalized 
non-Newtonian shear modulus. A low order polynomial 
ansatz for m suffices to model the feedback mechanism 
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of the cage-effect. We choose 

m(t) = T+W + W2$2 W) ' (3) 

Without shear, this model has been studied 
extensively^^. Increasing particle caging is mod- 
eled by increasing coupling parameters v\, v<i > 0, and 
the only effect of shearing is to cause a time dependent 
decay of the friction kernel m. The system loses memory 
because of shearing. The role of the transport coefficient 
(viscosity) 77 is played by the average relaxation time 
obtained from integrating the correlator. It is also taken 
to determine the stress: 

/>oo 

a = 7 V = 7 ( r ) = 7 / dt . (4) 
Jo 

At high shear rates, the memory function is strongly su- 
pressed, so that $ returns to a single exponential, and 
the high shear visosity of the model follows as rj^ = 1/T. 



D. Control parameter space and glass transition 
lines 

For the parameters of the model, the choice of generic 
values follows from previous considerations^* 2 ^. First, 
the parameter Y sets the time scale and determines the 
short time dynamics. The bare Peclet number for the 
model thus is given by Peo = 7/T. This parameter 
hence can be compared directly to the fast colloidal dy- 
namics determined by the radius of the particles and 
the short-time diffusion coefficient. Second, earlier stud- 
ies suggest to choose the two interaction parameters so 
that t>2 = «2 — 2 and v% = v\ + £/(^/v^ — 1), where 
vl = v$(y/4:/vZ - 1) w 0.828. Thus, the decisive pa- 
rameter, namely the effective volume fraction (j) of the 
particles enters the model only via e (</>). A glass transi- 
tion singularity lies at e = 0, where the long time limit 
<f>(t — > 00) = / jumps from zero for e < to a finite value 
/ > / c = 1 — 1/ yuf f° r e — anal 7 = 0. The parameter 
/ plays the role of the elastic constant Goo in this model. 



E. Flow curves of sheared suspensions 

The presence of a glassy arrested structure is equiv- 
alent to a frozen in part in the correlator or memory 
function; thus without shear $(i — > 00) = / > and 
m(t — > 00) = g > hold for e > 0. With shear a non- 
decaying part in m(t) is impossible, as 771(7:7 3> I) < 
+ v 2 )/(7£) 2 ; as a consequence, also $(f) always de- 
cays to zero. Memory is cut off at long times, and Eq. 
gives the most simple ansatz recovering this effect of 
shear advection in the microscopic equations 23 , and the 
obviously required symmetry in 7. The inset of Figure ^ 
shows the correlator for fluid and glassy states for systems 
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FIG. 1: Flow curves of the model, viz. dimensionless stress 
a versus 7/T, for two states that would be fluid respectively 
glassy without shear; e = ±£2 and e = ±£3 (where £3 = 
16 £2 = 0.0414). Straight lines with slope unity indicate the 
variation following the low shear 770 and high shear viscosity 

The inset shows the correlators of the schematic -model 
as function of rescaled time tT. The curve marked with a 
relaxation time r taken at ( I > (r) =0.1 corresponds to a fluid 
state without shear (e = —£2, 7 = 0). The curve marked 
by the long time plateau value / corresponds to a glass state 
without shear (e = +£2, 7 = 0). For increasing shear rates, 
log 10 (7/r) = -6, -4, -2, as labeled, the correlators decay 
more quickly at either value of £. 



at rest (7 = 0), and for sheared suspensions. For the lat- 
ter systems, the '-model predicts the speed up of the 
relaxation caused by increasing shear rates. Integrating 
over the correlators, as given in Eq. J3J|, leads to the vis- 
cosity which consequently exhibits shear thinning. The 
corresponding flow curves are shown in Fig. ^ Here the 
shear modulus a is plotted as the function of the reduced 
flow rate 7/T. Theory predicts an evolution from an (al- 
most) Newtonian fluid at weak coupling to a markedly 
non-Newtonian fluid at stronger coupling corresponding, 
in experiments, to higher volume fraction. This is seen 
from the characteristic S-shaped dependence of a on 7. 
At the glass transition, there is a discontinuous transi- 
tion from the dissipative fluid-like behavior to a yielding 
solid. A finite shear rate leads to a shear-melting of the 
glassy state. 

Fig. ^ presents the central result of theory. It gives 
the full scenario for the non-linear flow behavior of dense 
suspensions and relates it to the glass transition in these 
systems. Moreover, it predicts that flow curves obtained 
from glassy suspensions should present meaningful re- 
sults. This is due to the fact that the shear is expected to 
speed up the relaxation even at highest volume fractions. 
As an experimental consequence of this, no hysteresis is 
expected and flow curves present a well-defined probe of 
the dynamics of glassy systems. As a caveat, though, the 
condition needs to be recalled that the system is given 
enough time to reach the steady state, and that phase 
transitions and ordering phenomena are prevented. 
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III. COMPARISON OF MODEL 
CALCULATIONS AND EXPERIMENTAL DATA 
FROM MODEL SYSTEMS 



A. Thermosensitive latex particles 



A meaningful experimental study of the flow behavior 
of suspensions requires a system of particles that exhibit 
a rather small polydispersity and a high stability in the 
respective suspension medium. Moreover, the particles 
should interact in a purely repulsive fashion. This re- 
quirement is certainly given for the classical hard sphere 
suspensions used for the study of colloidal glasses so far. 
However, as an additional requirement, it should be pos- 
sible to prepare suspensions having volume fractions up 
to 0.63 in order to explore the region beyond the volume 
fraction of the glass transition. 

Recently, we showed that aqueous suspensions of ther- 
mosensitive latex particles meet these requirements^**^. 
The particles consist of a solid core of poly(styrene) of 
ca. 100 nm diameter onto which a shell of crosslinked 
poly(N-isopropylacrylamide) (PNIPA) chains is affixed. 
The particles are suspended in water and the PNIPA 
network in the shell is swollen at low temperatures (ca. 
10°C). Raising the temperature leads to an expulsion of 
this thermosensitive shell and the particles will shrink. 
The advantage of the thermosensitive suspension is obvi- 
ous: The effective volume fraction could be changed over 
a wide range by raising or lowering the temperature while 
keeping constant the weight fraction of the particles. In 
this way highly concentrated suspensions could be gen- 
erated in situ, that is, directly in the rheometer. No pre- 
shear was imposed on the system by handling or filling in 
a highly concentrated suspension. Moreover, any previ- 
ous history of the sample could easily be erased by lower- 
ing the effective volume fraction throuh raising the tem- 
perature. The thermosensitive suspensions hence present 
a novel model system by which volume fractions around 
and beyond the volume fraction of the glass transition 
become accessible without freezing-in non-equilibrium 
states caused by handling the suspension. Indeed, as 
shown in Refill 5 ., the shear viscosity rj of suspensions of 
these particles could be obtained over a large range of 
shear rates 7 in a well-defined manner. 

It should be noted that the thermosensitive particles 
present a well-studied system by now: The phase transi- 
tion within the microscopic network* 3 ^ was shown to be 
fully reversible 26-31 . Moreover, an analysis of these parti- 
cles by a combination of small-angle neutron and X-ray 
scattering revealed that the shell is well-defined and the 
particles exhibit a narrow size distribution 3 ^ 3 .. Their in- 
teraction in water is purely repulsive^! if the temperature 
is not raised over 30° C. 



B. Flow curves 

Theory states that flow curves of shear stress a ver- 
sus shear rate 7 refer to a well-defined stationary state. 
Close to vitrification, the parameters characterising the 
static structure can be subsumised into the separation 
parameter e. Increasing the interactions of the particles, 
brought about by increasing the volume fraction <fr then 
is described by increasing e from negative values in the 
fluid to zero, the point of the glass bifurcation, and be- 
yond, to positive values in the glass. 

Equation demonstrates that there is in principle 
only one parameter, the matching time to, that is re- 
quired to determine the flow curves and other steady 
state averages in the structural regime. This time scale 
contains all the effects of hydrodynamic interactions 
and other short-time phenomena not treated by theory. 
Hence, to needs to be adjusted by fitting the theoretical 
curves to experimental data. Let us stress again, that 
in principle one time to suffices to describe all different 
experimental measurements on a sheared dispersion, for 
all volume fractions close to its glass transition, and for 
small shear rates. 

An 'idealized' analysis of a complete experimental flow 
curve would thus proceed via 

C -> 7 Voo + ^struct. (£, 7*o) 

~ ^struct, (e, 7*o ) f° r |7*o| < 1, N < 1 , (5) 

where r]^ describes the flow curve at high shear rates. 
At small shear rates, the two fit parameters 7]^ and to 
remain as unknowns for the description of the structural 
region in the flow curve, as cr s t xuc t.(e, x) is determined by 
the static structure at the transition (viz. the critical val- 
ues v? of the vertices in Eq. ©). Both parameters, 7700 
and to, are influenced by the physics at high shear rates 
and short times, that is dominated by hydrodynamic in- 
teractions. 

A representation where the subtleties of the non-linear 
flow curves are revealed most clearly is given when plot- 
ting the stress versus shear-rate. Figure [21 shows the 
flow curves of the thermosensitive particles for interac- 
tion strengths close to their glass transition. The data 
have been taken from Refill 5 .. Decreasing the tempera- 
ture T swells the particles such that the effective packing 
fraction can be thought to increase. A region for low 7, 
where a depends strongly on density, can be seen apart 
from one at high 7, where a smaller variation is found. 
It is the region at low shear rates which is treated by 
the theory, i.e. where structural dynamics dominates the 
flow curves. Theory suggests to plot the stress a ver- 
sus shear rate 7, instead of viscosity 77 versus 7. This 
provides a direct comparison with Fig. The lack of 
straight pieces in Fig. [21 indicates the absence of true 
power-law shear-thinning, r\ cx j~ x , which would show 
up as a oc j 1 ~ x . 



5 




-2-10 12 3 



iog 10 f/s" 

FIG. 2: Flow curves, viz. stationary transverse stress a as 
function of shear rate -y, for the thermosensitive latices close 
to glassy arrest at various temperatures T from top to bottom 
as denoted. The temperatures may be thought to correspond 
to different effective packing fractions. 

C. Elimination of effects from hydrodynamic 
interactions 

The ideal analysis of the structural part of experimen- 
tal flow curves, which only requires to to be matched at 
the transition point, is hindered by (i) inevitable quan- 
titative errors of the theory in calculating various con- 
stants, like e.g. the critical packing fraction, the yield 
stress at the critical point, the transversal elastic con- 
stant, and other quantities that are determined by the 
dispersion structure. Obviously, these quantities can- 
not be calculated within schematic models where a small 
number of vertices {v{\ replaces the structural informa- 
tion. In a schematic model analysis, therefore, the overall 
stress amplitude has to be fitted; 

Struct. = (T * a thC °- ({Vi}, 7*0 ) . (6) 

The time scale to is easily determined within the 
schematic model, and may then for convenience be elim- 
inated in favor of the intrinsic decay rate T of the model. 
We adopt this convention, as T can be read off more eas- 
ily from the (7(7) curves. 

The ideal analysis must also require that (m) param- 
eters like (To, ?7oo, and to that are not treated by the 
model are constants. In particular, 77,30, and to are re- 
lated to hydrodynamic interactions. To accomodate for 
these unknown parameters and their non-negligible den- 
sity dependences, an analysis using schematic models can 
proceed via relaxing the restriction that all parameters 
except for e are constant. The data shown in Fig. [3] 
exhibit a density dependence of the stress at high 7 
which is not contained in the i^j '-model; comparison 
with Fig. n snows that the model leads to a constant, 
(Too = 7?/oo = 007/T, which thus does not constitute a 
full description of the hydrodynamic effects. Thus in Eq. 
©, a temperature dependence needs to be included in 



the parameter <Too(T). As discussed in Ref»2&, we expect 
that this needs to be accounted for by including a tem- 
perature dependence of (Jq. In order to keep the number 
of temperature dependent parameters as small as possi- 
ble, in the following analysis all temperature dependent 
corrections beyond the model are assumed to arise from 
00 (T), the temperature dependence of the overall stress 
prefactor. The final expression used for the data analysis 
with the -F^p-model thus becomes 

<j = a (T)*a thc °-(e(T),j/Y) (7) 

where, e captures the temperature dependence of all ver- 
tex parameters, as has been discussed repeatedly in the 
literature for this model without shear. 
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FIG. 3: Stress data from Fig. [5] rescaled by x — ao(T = 
10°C)/a (T) to agree at high shear rate in order to eliminate 
corrections from i.a. hydrodynamic interactions; rescaling 
values x from top to bottom as indicated in the legend. Solid 
lines give fits with the schematic model to the flow curves at 
small 7 where the structural dynamics dominates. The pa- 
rameters are a (T = 10° C) = 4.69 Pa, T = 280s" 1 , for all 
curves, and separation parameters e(T) as given in Table I. 
The glass transition temperature lies close to T c — 11 C. 



D. Data analysis 

The parameters obtained from fits of the schematic 
model curves to the experimental data are given in table 
I. The fits are shown in Fig. [3] Adjusting the parameters 
starts with the data sets closest to the glass transition at 
e = 0, where the overall scales <To and Y can be found 
as the flow curves show the strongest variations in cur- 
vatures, and then proceeds to the temperatures farther 
away. 

The transition to a yielding solid at e(T c ) = is found 
to lie around T c as 11°C, where the critical value of the 
yield stress is er+ = 0.44 Pa. It needs to be noted that the 
transition is not brought about directly by the change of 
temperature but by the change of the volume fraction of 
the particles with temperature. For lower temperatures, 
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T in°C 


10 


15 


20 


25 


30 


e 


0.010 


-0.037 


-0.108 


-0.20 


-0.315 


(To in Pa 


4.69 


2.93 


2.35 


1.62 


1.09 


?7oo in 10 _2 Pa s 


1.67 


1.04 


0.84 


0.58 


0.39 


770 in 10~ 2 Pa s 


00 


341 


27 


5.7 


1.7 



TABLE I: Parameters, e(T) and <x (T), for the fits of the flow 
curves in Fig. fusing the -model. The viscosities 770 and 
7700 are calculated from the fits using P = 280s -1 . 



the yield stress increases quickly, a + (T — 10°C) = 1.3 
Pa. At the fluid side of the transition, the zero shear rate 
viscosity 770 increases by more than two decades, while 
the high shear rate viscosity rj^ changes by less than a 
factor 5. 77^ is not predicted by the present model, it 
arises from hydrodynamic interactions, but it is taken 
into account by by varying <jq(T) with temperature. 




IV. CONCLUSION 



A quantitative comparison of a first-principles ap- 
proach to the nonlinear rheology of dense colloidal 



system; 



23.24 



with experimental data obtained in model 
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dispersions^ has been given. The employed schematic 
F 12 -model of the nonlinear rheology has been obtained 
after simplifying steps based on the microscopic mode 
coupling theory. The speed up of the structural relax- 
ation brought about by shearing the suspension is the 
central mechanism considered which causes shear thin- 
ning and yielding behaviors, ft enters the schematic 
model via a time-dependent suppression of long-term 
memory. 

The comparison with the experimental data demontrates 
that this model captures the essential features of the flow 
behavior of concentrated suspensions: i) the strong shear 
thinning with increasing shear rate, and ii) at low shear 
rates the transition from a Newtonian liquid to a soft 
yielding solid. Using the usual assumption that the hy- 
drodynamic interactions can be described solely in terms 
of a high-shear viscosity 77,30 a full description of the ex- 
perimental data by Eq. has become possible (see Fig. 
|3J. The comparison demonstrates that the strong raise 
of the viscosity with increasing volume fraction can be 
fully explained by the structural arrest of the particles 
when approaching the volume fraction of the glass tran- 
sition. It is hence evident that mode coupling theory 
that provides an excellent description of the dynamics of 
quiescent suspensions yields also a quantitative explana- 
tion of flow curves observed for suspensions subjected to 
a steady shear field. 



FIG. 4: Viscosity data and corresponding schematic model 
fits taken from Fig. |3] The second Newtonian plateau in the 
fit curves increases with lowering temperature because of the 
rescaling factors x. Curves from top to bottom as labeled. 

Replotting the data and model results in the classi- 
cal flow curves showing viscosity versus shear rate. This 
enlarges the variation along the ordinate and thus sup- 
presses the slight deviations of the fits from the data. 
Figure 0] shows the curves without rescaling so that the 
temperature dependence of the high shear rate viscosity 
7700 becomes apparent. 
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